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Abstract 

This paper is an electronic application to my set of lectures, sub- 
ject:'Formal methods in solving differential equations and constructing models 
of physical phenomena'. Addressed, mainly: postgraduates and related read- 
ers. Content: a very detailed discussion of the simple model of interaction 
based on the equation array: 



d 2 q(t) 
dt 2 

d 2 u(t,x) 2 d 2 u(t,x) 



(«(*) - Q(t)) + Mt) 

4rfc8(x-xo)(Q(t)-q(ti) + fi(t,x) 



dt 2 dx 2 
Q(t) = u(t,xo) 



Besides, less detailed discussion of related models. Central mathematical 
points: dAlembert-Kirchhoff-like formulae. Central physical points: phe- 
nomena of Radiation Reaction, Braking Radiation and Resonance. 
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Introduction. 

A Harmonic Oscillator Coupled to an One-Dimensional Scalar 
Field. 

In this paper I will discuss several models of an one-dimensional harmonic oscillator 
coupled to an one-dimensional scalar field. Primarily I am interested in the model 
described by the equation array u 

d 2 q(t) 



dt 2 

d 2 u(t,x) 2 d 2 u(t,x 



c: 



n 2 (q{t)-Q{t))+f {t) 

4 7 c6{x - x )(Q(t) +h(t,x) 



dt 2 dx 2 
Q(t) = u(t,x ) 



One can say that this is a model of a point interaction. The considered model is not 
standard, if one means classical or quantum field models. In that field, a standard 
looks rather like this: 



d 2 q (t) 



= -n 2 q (t)+j 1 Q 4 ,(t) + f (t) 



dt 2 

d 2 (j){t,x) _ 2 d 2 4>(t,x 



c 



dt 2 dx 2 
Qt(t) = (t>(t,x ) 



+ 4~f 2 c5(x - x Q )q (t) + fi(t,x) 



After indicating d'Alembert-Kirchhoff-like formulae for solutions to these 
systems I briefly compare them and discuss phenomena of Radiation Reaction, 
Braking Radiation and Resonance. n 

After that I will discuss an abstract analogue of the former system: □ 

q = -n 2 (q-Q) + f (t) 

u = Bu- 4rf e (s a ,t^ ■ (Q - q) + fi(t) 

Q = Q{t) =< l\u(t) > 



where I use a P.A.M. Dirac's "bra-ket" syntax and suppose that q and Q are usual 
(one-dimensional) functions of t : 

q = q(t) , Q = Q(t) , 

B is an abstract linear operator, I is a linear functional, {u(t)} t and {5 a j.x }t are 
families of abstract elements; of course the type of S a , t ,x must be the same as one 
of u(t). 

One can rewrite the above equations as follows: 

q\ = ( n 2 <l\ 2 \(q \ ( f Q (t) 

U ) \ ^cSa,t,x < l|i B-Aj c S a:ttX0 <l\ 2 ) \ u ) \ A(t) 

Qo\_f 71 < l\2 \f qo \ ( fa{t) 

^2,oS a ,t, X0 < l|i B J\ 4> J \fi(t) 



1 47cp = k , Q 2 = k/M ; the constants of the model mean, e.g.: c = propagating waves velocity, 
p = 'a density' of the field, k = elasticity constant, M = mass of the particle. Of course, we 
assume c > . 

2 the proper modifications for the latter system will be evident 
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where the subscribts 1 and 2 in < and < • ■ ■ I2 mean that arguments of 

< • • • |i are elements of the first component of the vector f ^ ~\ resp. ( ^ ^ and 

arguments of < ■ ■ ■ I2 are elements of the second component of the suitable vector. 
Of course, 

< l\q > = < >=q, < l\q >=< l\iq >= <7o 

Normally I will suppose that 8 a ^.x a does not depend on t, i.e. is constant in t. 
In that case I will write 5 a ,x instead of 5 a ,t,x - 

In the next paper I will suppose that 5 a ^.x lS constant in t. 

The subject will primarily be resolvents formulae i.e. the formulae that re- 
solve the equation array: 

u ) \ 4rf c 6 a , xo < 1 1 x B-A-f c 6 a . Xo < l\ 2 J \ u J \ w 2 

resp. 

qo \ _( 71 < l\ 2 \ f qo \ = ( Wl 

4> J V ^I2,c5 a ,x < l|l B ) \ (j) ) \w 2 

Of course, here all quantities, q,u,wi,w 2 ,S a . XQ ,l, etc. are supposed to be con- 
stant in t. n 

I will discuss Donoghue-Friedrichs-like solutions to the system til and resp. the 
phenomenon of Resonance and notion of the Second Sheet. 



3 i.e., I will discuss resolvents formulae of the system 
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1 Models of a Point Interaction of an only Oscil- 
lator with an only one-dimensional Scalar Field 

In this section we fix measure units and let x be dimensionless position parameter, 
i.e., 

physical position coordinate = [ length unit ] x x + const . 
Otherwise a confusion can ocurr, in relating to the definition 



J —i 



S(x - x )f(x)dx = f(x ) 



We assume the standard foramalism, where 

dl+(x-x ) 
5(x-x ) = 

and where 1 + stands for a unit step function (Heaviside function): 

1 , if ^ > 0, 



o , if e < o , 



1.1 First Model. D'Alembert-Kirchhoff-like formulae 

Recall that standard D'Alembert-Kirchhoff formulae read: if 

d 2 u id 2 u , . . , ,. . , . 

-Qi2= C Q^2 + f' u = u(t,x), f = f(t,x), (*) 



and given initial data, u(s, •) and 9u ^'& ) 



then 



U = u(t,x) = ^J(f(T,X + c(t-T))-f(T,X-c(t-T)f)dT 

+ u (t,x) 

Uo{t,x) = c + (x + ct) + C-(x — ct) 

= ^(u(s,x + ct) + u(s,x — ct)^j 

■ (u(s, x + ct) — ii(s, x — ct)j 



and where / , u stand for any functions defined by 

df(t,x) . dZ( s ,t) fdu(t,zy 



~ n ' h dt ~\ dt J 



dx J v ' " d£ 

Note that 

/(r, x + c(t — t)) — /(t, x — c(t — t)) , ii(s, x + ct) — ii(s, x — ct) 

do not depend on what the primitives are which one has chosen!!! Moreover, we 
need only / , u|t =s and not / , ?i| t=s themselves! 

Recall that standard relation for an one-dimensional harmonic oscillator under 
an external force f (t) is this: 



d 2 q(t) 

dt 2 



&q(t) + f (t) 
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In this section I will firstly discuss a system described by 
8 2 q(t) 



dt 2 



-n 2 [q(t)-Q(t))+f (t) 



d 2 u(t,x) _ ^d 2 u{t,x 



dt 2 dx 2 
Q(t) = u(t,x ) 



4jcS(x - x )(Q(t) - +fi(t,x) 



It means in patricular that I will take 

/ = -A 1C 8(x - x ) (g(t) - <?(*)) + fi(t, x) 
For such an / I conclude that 

f = -A 1C l+{x - x )(Q{t) - q{tj) +fi(t,x), 



and then I infer that 



u(t,x) = -1l j (l+(x + c(t-T) -.t ) - \+{x-c{t-T) -x )^(q(t) -g(T))dr 

+ Yc j s X + C ^ ~ T ^ ~ x ~ c ^~ T ))) dr + u o(*, s, x) 

u(t,x) = _2 7y (l+(x + c(t - t) - x ) - l+{x - c(t - t) - x fj(u(T,x ) - q(rfjdT 

+ ^c/ (fl( T > X + C {t- T )) ~ fl(T,X-c(t-T))^dT + U (t,S,x) 

Denote now, to be more concise, 

u 01 (t,s,x) := Y c \ a (/i( T ' a: + c ( <_T )) ~ fi(T,x - c(t- T))^dT + u (t,s,x) 
and then rewrite the recent relation as following: 

u{t,x) = -27 J (l + (x + c(t-T)-xo)-l+(x-c(t-T)-xofj(u(T,xo)-q(TfjdT 

+ u Q1 (t,s,x) 
I have now seen: given q and uqi, then, 



in order to obtain u(t, x) I need to obtain ONLY u(t, xo) 



After this observation use the last formula for u(t, x) and then obtain 

u{t,x ) = -2 7 ^ (l+(c(t-T))-l+(c(T-t)))(«(T,a:o)-«(T))dT 

+ U l(t,S,X ) 

Therefore, because 

l+(c(t - r)) = 1 , l+(c(r - i)) = for t > t 

one can obtain 
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i.e., 



Write it as 



u(t,x Q ) = -2 7 (u(t,x ) - q(r)jdT + u i(t,s,x ) 
Q(t) = -2 7 ^ (Q{T)-q{T)yT + u 01 {t,s,x ) 

Q(t) + 2~/ [ Q{r)dT = 2 7 / q(T)dT + u i(t,s,x ) 

J S J S 



Now recall that 

d 2 q{t) 

at 2 

and for the moment denote 



-Q 2 [q(t)-Q(t))+f (t) 



Qo(t) := u (t, s,x ) = c + (x + ct) + C-(x - ct) 

Qoi(t) ■= Uoi(t,s,x ) = Yc\ s (/ 1 ( T ' a;o + c ( t ~ T )) ~h{T,x -c(t-T)fjdT + Q a (t) 
Then obtain 

Q(t) + 2j [ Q(r)dT = 2 7 / q(r)dT + Q Q1 (t) 

J s J s 

= -n 2 (q-Q) + f (t) 

We have now obtained an equation array for Q and q. The next step is to find 
an insulated equation for Q and one for q. We begin to search for the equation for 
q. For this purpose, we apply, at first, the operator 



( I + 2 lj s - dT ) 



to the latter equation, i.e. to the equation 

(Pq 

dt 2 

Then we infer 

32„ rt e>2 



i)Iq = -n 2 (q-Q) + f Q (t) 



+ 2 (q + 2 7 ^ Q(r)dr) 
+ (/o + 2 7 J /o(r)dr) 

+ tt 2 (2 1 q(r)dr + Q 01 (t)) 
+ (/o + 27^ /o(r)dr) 



2„ /•« fl2 



at 2 7 



£^$r dT = -n 2 (q + Qoi(t)) + (.fo + 2jJjo(T)dr) 
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The most recent equation for q rewrite as 



d 2 Q „ 



j ^rdT-n 2 q + n 2 Qoi(t) + f (t) + 2 1 J^ f ( T )d T 



We can stop at this equation, or, observing that 



t d 2 q(r) dT _ dq(t) dq(t) 



<9t 2 



dt dt 



we can stop at that: 



dt 2 =- 2 ^- n ^ + ^dt 



+ Q 2 Qoi(*) + /o(*) + 2 7 / / ( r )dT 



/' 



Some people prefer to write such an equation as following: 



^ + 2 7 ^ + tf)q(t) = 2 7 ^|^ + n 2 Q 01 (t) + fo(t) + IjJ* h{r)di 



Remark, whatever equation we take, we see an interesting detail: the equation is 
not ordinary differential equation, if we follow standard terminology. The case 
is because of the term 

dq(t) 



2 7 - 



dt 



In the next subsections, we will return to this factor. However, the machinery of 
the ordinary differential equations does here quite for. For the time being, we turn 
to describing Q(t) and u(t, x). 

After q(t) is found, we can determine Q(t) , at least formally, by solving 
Q(t) = -2 7 ^ (q(t) - q{T)y T + Q 01 (t) 



or 



dQ(t) 
dt 



= -2 7 (Q(r) - q(r) 



dQoijt) 
dt 



Q(s) = Q 01 (a) 



Thus we have already reduced our model, a model of an oscillator coupled the 
a scalar field, to a pair of linear 'ordinary' differential equations. Nevertheless we 
want to continue to analyse the matter and we now go seeking another relationships, 
which would simplify calculations of q, Q and u. 

At first, we will obtain another insulated equation for Q, differently and in a 
different form. 

We have 



Q(t) + 2 7 f Q(r)dT 

J s 



2 7 / q( T )dT + Q i(t) 



^| + o 2 ? = n 2 Q + f (t) 
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i.e., 



(Q(*)-Qoi(*)) +27^ (Q(r)-Qoi(T))dr = 2 7 ^ ( 9 (r) - Q i(r))dr 

(^ + r! 2 ) (J = n 2 Q + / (t) 



Then we infer 
9 2 



+ f! 2 ) ((Q(i) - Qoi(t)) + 2 7 / (q(t) - Q i(r))dr 



(S +fi2 )( 27 /* ( 9(T) ~° 01 (T) ) dT ) 



and 



(^2 + fi2 ) ((<?(*) - Ooi(t)) + 2 7 ^ (g(r) - g i(r))dr) 

= {w + n2 )i 21 f s q{T)dT ) -{^ + n2 )i 21 f s Qm{T)dT ) 



Note 
Then 

32 



W2 l q{ r) d r = - q{ t)=J s — q{T ) d r + - q (t) 



t=s 



(^2 +" 2 )((0(*)-Qoi(t)) + 2 7 ^ (Q(r)-Q i(r))dr) 



Then 

a 2 



(^2 +" 2 )((Q(*)-QoiW) +27 y (Q(r)-Q i(r))c 

= 2 7 jT (fi 2 Q(r) + /o ( r ))dr + 2 7 J^(t) | ^ - (7^ + ^) (27^' Qoi(r)dr) 



Then 



+ ° 2 ) ((<?(*) - Qoi(*)) + 27^ (<3(r) - Q m {T))dr 



= 2 7 ^ (fi 2 Q(r)-r! 2 Q 01 (T))dT + 2 7 ^^) 
+ 2 7 / / (r)dr 



t=s at 2 



(2 7 ^ Qoi(r)dr) 



Then 



(S + (° w ~ ° oiW ) + (S + " 2 ) 27 / t ( Q(r) ~ Qoi ^) dT 
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2 7 ^ (fi 2 0(T)-^ 2 Qoi(T))dT + 2 7 ^g(i)| t=s -^(27^ Qoi 



(r)dr 



+ 2 7 



I h{r)di 

J s 



Then 

(&_ 

\dt 2 



tf)((Q(t) - Qoi W) + ^2 7 / (0(r) - Qoi(r))dT 



= 2 7 - ? (t) 



t =. at 2 



(2 7 ^ Qoi(r)dr) +2 7 ^ / (r)d 7 



Then, finally, 



(£ + 2 4 + fi2 )( Q(<) -°° l(<) ) 



2 7 



0g(t) 



at 



27 — ^ — + 2 7 J fo(r)di 



On the surface, this equation appears to be a second order ordinary differential 
equation. It is not exactly the case. We may not arbitrary take the initial data for 
Q(t). The proper ones are these: 



(Q(s) - Qoi(s)) =0. 



d(Q(t)-Q 01 (t) 



dt 



101 



(s) 



Thus, we have already obtained two simple 'ordinary' differential equations for 

q(t) , Q(t) 

Now let us analyse the expression 

u = u(t,x) = -2 7 J (l+(x + c(t - t) - xq) - l+(x - c(t - t) - x q )^(q{t) - q(r)jdT 
+ u i(t,s,x) 

We have: if ct ^ ct + (x — x ) and ct ^ ct — (x — x ) then 

l + (x + C(t - T) - X ) - l + (x - C(t - t) - Xo) 

_ J 1 , CT < Ct + (x — Xq) / 1 , Ct — (x — Xq) < T 

I 



Ct + (x — Xq) < CT 



{ 



CT < Ct — (x — Xq) 







CT < Ct 


-\x 


- x \ 




■1 


Li 


ct — \x - 


- x 


< CT < ct + \x 


- %o\ 






ct + \x - 


- X 


< CT 








T <t- 


\x- 


x \/c 




■ 1 




t-\x- 


xo\l 


C<T <t+\x- 


- X \/C 






t + \x - 


Ml 


C<T 





Hence 
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u(t, x) 



/t-\x-x \/c 
\Q{t) - q(r))d,T , if s < t - \x - x \/c 



-2 7 



rt+\x-xa\/c 



if t — \x — xo\/c < s < t 
iit<s<t+\x — xq\/c 



[Q(T)-q{T)JdT , if t+\x-x \/c< s 

+ u 01 (t,s,x) 



Besides, remember that 



Q(t) + 2j f Q(t )dr = 2 7 / q{r)dT + Q i(t) 

J s J s 

If we take this factor into account, then we deduce, finally, 



u(t, x) 



Q(t - \x - x \/c) - Qoi(t - \x - x \/c) , if s < t - \x - x \/c 

, if t — \x — xq\/c < s < t 

, if t < s < t + \x - xo\/c 

-Q{t + \x - xq\/c) + Qoi{t + \x - xq\/c) , if t + \x - x \/c < s 

+ u i(t,s,x) 

Now substitute anywhere u (t,s,x ) for Qo(t), and u i(t,s,x ) for Qoi(*) ; an d 
resume. 
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1.1.1 Putting It Together 

We have discussed the system 
d 2 q(t) 



dt 2 



= -fi 2 (g(t)-Q(t)) +/ (t) 



d u{t,x) 2 d u(t,x) 

dt 2 dx 2 
Q(t) = u(t,x ) 



- AjcS{x - x ) (0(t) - + fi(t, x) 



and concluded that: 



(w + 2l ¥t +n2 ) q(t) = ^^\ t=s + n 2 uoi(t,s,xo) + Mt) + 2 1 j\fo(T)dT 
Q{t) = -2 7 ^ (q(t) -q(r)yT + uoi(t,s,xo) 



An insulated equation for Q(t) is this: 



= 2 7 - 



0uoi(t,s,a:o) 



0t 

The proper initial relations are these: 



t=s Ot 



t fo(r)dr 

J s 



. , d(Q(t) -uoi{t,s,x )) , , 

(Q(s) - u 01 (s, s, x )J = , — ^ — '- ^ = 2j[q(s) - u 01 (s, s, x )J 

For u — u(t, x), we have concluded: 
ft 

u(t,x) = -2 7 y (l + (x + c(t-T)-x )-l + (x-c(t-T)-x )^(Q(T)-q(T)^dT 
+ uoi(t, s, x) 
Other expressions for u = u(t,x) are these: 



u(t,x) = 



Q(t-\x-x \/c)-u 01 (t-\x-x \/c,s,x ) , if s < t - \x - x \/c 

, if t - \x - xq\/c < s < t 

, if t < s < t + \x - xq\/c 

-Q(t + \x - x |/c) + uoi(t + \x - xo\/c, s,x ) , if t + \x - x \/c < s 



+ u i(t,s,x) 
where uoi is defined by the relations 

Uo{t, s,x) := c + (x + ct) + c-(x — ct) 

= -( u (s,x + c(t-s))+u(s,x-c(t- S ))) + -j^^(^^) <li, 
u 01 (t,s,x) := Y c \ a (M T > X + C (t- T ))- h( T > x - c ( t - T ))) dT + u o(t,s,x) 



t=s 
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In the following subsections, we will mostly discuss the case where 

c=l,s = 0,a;o = 0,/i=0. 

it will be convenient to have rewritten some of the recent formulae in the proper 
way. 

(S + 2 ^Ft + n2 ) q{t) = 27 ^f? L + n2uo{t > °' o) + m + 27 Si fo{T)dT 



Q(t) = -2 7 ^ t (Q(T)- <Z (r))dr + Mo (t,0,0) 



An insulated equation for Q(t) is this: 



(^ + 27 ^ + " 2 )( Q(<)_Uo(i '°' 0) ) = 27 







dq(t) 



dt 



t=o 27 01 



du Q (t, 0,0) 



+ 27 / fo 
Jo 



The proper initial conditions are these: 



, N d(Q(t)-uo(t,0,0)) 

(Q(0) - «o(0, 0, 0)) = , ^ t=Q = 2 7 (g(0) - uo(0, 0, 0)J 



Expressions for u = u(t, x), are these: 
u(t, x) = 



Q(t- \x\) -u (t- \x\,0,0) , if0<t-|a;| 

, ifi-|a;|<0<t 

, ifi<0<i+|a;| 

-Q(t+ \x\) +u (t+ \x\,0,0) , ift+|a:|<0 



where 



+ u (t,0,x) 



Uo(t,0,x) := c + (x + t) + C-(x — t) 



1 / \ 1 f x+t 

-[u(0,x + t)+u(0,x-t))+- J ^ ( 



du(t,£). 
dt 



t=o 
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1.2 A Standard Model. D'Alembert-Kirchhoff-like formulae 

The considered in the previous subsection model is not standard, if one means 
classical or quantum field models. In that field, a standard looks rather like this: 

d 2 q (t) 



dt 2 

d 2 4>{t,x) 



-n 2 q (t)+ 11 Q 4> {t) + f (t) 



,d 2 cj>(t,x) 



dt 2 dx 2 
Q<l,{t) = <f>(t,x ) 



+ 4-y 2 cS(x - x )q (t) + fi(t,x) 



In the latter situation one can obtain that 

rt 



Q<t>{t) = 272 / qo(r)dT + 4>ai(t, s, xq) 

J s 

for the proper <poi(t, x) . Next, with the reasons of the previous subsections: 




and 



4>o\{t, s, x )j 

d 2 



2 l2 (^ + n 2 ) j\ (r)dr 

272 / (jJ +n2 ) 9o(r)rfT + 272 -| 9(<) 
272 f {iMt) + /o(r))rfr + 2 7 2 J^(i) 

27271^ (Q<p( T ) - (t>Ol(T,S,X )^dT 



+ 272 J (li<t>oi(T,s,x ) + fo(r)^dT + 2"f 2 ^:q{t) 



'dt 



Thus 



(g^ + n2 ) (Q*(*)-0Ol(t,«, 


&o)) 






= 27271 




(t) - (t>oi{T,s,x )^dT 




+ 2 72 


S, ( 7i< 


\ d 
Poi{t,s,x ) + fo(T)jdT + 2 72 — q(t) 


t=s 



We can reduce these equations to ordinary differential equations, but these 
ones become third order equations, 



d 3 . n2 d \ d<p 01 (t,s,x ) df (t) 

di~ 21211 ) qo{t) = 71 at + — 



( 



(^3 + ^ 2 ^ + 27271) 0oi (t,s, a*)) = 272 (71^01 *, &o) + /o(t)) 
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subject to the supplementary initial relations 

d 2 



G* +na )*< 4 >L = /o(s) 



(^2 +" 2 ) 001 (*.*,*(>)) 



2^q{t) 



As for (j)(t,x), one can obtain that 

</>(t,a;) = 2 72^ ^l+(a; + c(t - t) - a?o) - l+(a; - c(i - t) - aro)^go(T")dT 



Other expressions for </>(i, x) will be these: 

rt— \x — xq\/c 



272 / 
J s 



<j)(t,x) 



qo(r)dT , if s < t — \x — xq\/c 



rt+\x — xo\/c 



f 

-272 / qo{r)dT 
+ <f>oi(t,s,x) 



if t — \x — xq\/c < s < t 
if t < s < t + \x - x \/c 

if t + \x — xq\/c < s 



<P(t,x) = < 



Q^(t-\x-x \/c)-<j>oi(t-\x-XQ\/c,s,xo) , if s < t - \x - x \/c 

, if t — \x — xq\/c < s < t 

, if t < s < t — \x — xq\/c 

-Q < t,{t+\x-xo\/c) + 4>Qi{t+\x-x \/c,s,x ) , if t + \x - x \/c < s 

+ <f>oi(t,s,x) 

where <poi is defined of course in the same manner as in the previous subsection, by 
the relations: 



b (t,s,x) := C0 i+ (x + ct) + c$-(x - ct) 

= \ x + c(t- a)) + <f>(a, x-c(t- «))) + ^ £ 

01 (t,s,x) := ^ J (ji(T,x + c(t - t)) - fi(r,x - c(t - r))^dT + (j> {t,s,x) 



1 f x+cit - s) ^ d(j)(t,0 

x-c(t-a) & 
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The case where 
is this: 



c=l,s = 0,a;o = 0,/i=0. 








= 2 72 7i J ( 


Q (r) -Mr, 0,0)) dr 


+ 2 7 2 / 

Jo 


(yiMr, 0, 0) + h(r))dT + 2 7 2^(i) f _ Q 



Reducing to ordinary differential equations: 
(^3 +" 2 ^ + 2 7 27i)(Q0W-0o(i,O,O)) = 



.dt 3 ' " dt 
supplementary initial relations: 



afo(t,o,o) , df Q (t) 

71 — — + — 

2 72 ( 7 i0o(t,O,O) + / o (i)) 



(& + ° 2 HC = /o(0) 



(^2 +^ 2 )(Q4*)-^o(t,o,o) 



t=0 



2 72 |,W 



t=0 



As for <f>(t, x) 



<j)(t,x) = < 



Q*(t-M) -&>(*- N.o,o) 
o 
o 



+ 0o(*,O,ar) 
where 0o is defined by the relations: 

0o(*,O,a;) := C0 i+ (a; + i) + c$-(x - t) 



if < t- \x\ 
if f - \x\ < < t 
if t < < t + \x\ 



-Q4,{t+\x\) + (!>o(t+\x\,0,0) , ift+|a;|<0 



1 / \ 1 f x+t 

-(<l>(0,x + t) + <f>(0,x-t)j+- J ^ ( 



dt 



t=o 
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1.3 Particular Cases. 1. Radiation Reaction, Braking Ra- 
diation 

We are all familiar with the fact that any solution to any linear inhomogeneous 
equation, whatever its nature, is a sum of a solution to the associated homogeneous 
equation plus arbitrarily taken and fixed solution to the former linear inhomoge- 
neous equation, isn't it? Just let us now consider these two cases separately, the 
case where equations are homogeneous and then inhomogeneous. 

We start out emphasising that the homogeneous equation array, connected to 



(w + 2 4 + n2 ) q{t) = 2l ^rl=o + n2uo{t ' 0,0) + fo{t) + 27 / " fo{T)dT 



dq(t) 



( d 2 „ d 
I 1- 2 7 — 



\dt 2 
is exactly 



dt 



fl- 



at 

Q{t) - u (t,0,0) 



= 2 7 



dq(t) 



dt 



27 



du Q (t, 0,0) 
dt 



27 / k{r)dT 
Jo 



\dt 2 



(w +2 4 +n2 ) qit) = 27 

+ 2 7 ^ + fi 2 )(Q(t)- Wo (i,0,0)) = 2 7 



dq(t) 



dt 
dq(t) 



dt 



t=0 



and NOT 



\dt 2 



+ 2 7 ^ + fi 2 )(g(t)- Uo (i,0,0)) = 



The difference is a rank one term 2 7 



dt 



t=0 



This detail allows us to apply 



usual machinery of finite rank perturbations theory. Thus, having put 



)2 _ ^2 



- .- \/n 2 

and having taken into account the supplementary initial relations: 

dQ(t) 



0(0) = 0, 

one can show that 
q(t) = e-^ (cos(Q 7 t) + ~- 



dt 



t=o 



(?(o) 



2 7 g(0) 



27 dq(t) 



Q 2 dt 



t=o 



_ 7t sin(Qyt) dq(t) 



at 



t=0 



n 2 dt 



Q(t) 



- «-7t 



(cos(fljt) 
2 7 dq{t) 



sm(f2 7 i) \ / 2 7 dq(t) 



n 2 dt 



t=o 



n 2 dt 



_ 7t «nflV) 27g(0) 



4 This term, as a function of t , is a fixed function of t, e.g. 1, multiplied by a CONSTANT 
depended on q, i.e., by a fixed functional of q. Using the Dirac's syntax, this term can be written 

as a >< b\ with a >= 1 and < b\q >= 27 
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and then 

u(t, x) 



u(t, x) 



= e 



-7(t-| 



X D (cosn 7 (t - \x\)+ 7 



sinCt 7 (t — \x\) 



ft 



)( 



2 7 0g(t) 



ft 2 St 



J 



+ e _ 7(t _ |x | )g m0 7 ft-N) 27g(Q) 



2 7 0g(t) 



+ at 

, if o > t 



if < t 



t=o 



Notice, 



g(t)-Q(t) = e -^ COS (^(0) + e ^ tf ^(^[ =o -2 7g (0)) 



and then 



q(t) - Q(t) ->■ as i ->■ +oo , 7 > 



If we concentrate now on the system's behaviour at large t , a mathematical 
detail calls attention. We observe: 



q(t) -> 2 7 



S«(t) 



St 



as t 



+00 . 



t=o 



Nevertheless, the limit function 



<7oc(t) = 27 



dq(t) 



dt 



t=0 



is NO solution to 



every time that 



d 2 q(t) 

dt 2 



-2 7 -^ - r» 2 g (t) + 27. 



27 



9«(t) 



at 



t=0 



dt 



because the q = q^ is a constant, hence its derivative is zero but we have need of 



dq(t) 



dt 



dq 



_d_( dq{t) \_ 
~ dtV 1 dt t=o) ~ U 



t=o 



dt 



Similar phenomena, one can detect them in the electrodynamics of moving charges. 

Now, we are directing our attention to the fact that the oscillator moves with 
damped amplitude: 

\q(t) — qoo I < e~ 7 * • const . 

The reason is very plain. When the particle (oscillator) begins to move, it makes 
some region of the held 'move': the oscillator generates, emits waves, -one is used 
to saying: 'the oscillator radiates'. 

Now then, the oscillator radiates. It entails some energy expenses and thus 
the field, made move, eventually makes the conduct of the oscillator change. This 
phenomenon is said to be a radiation reaction. 

Next, since the oscillator has emitted waves, waves run away and brings away 
some portions of the oscillator energy, braking the moving of the oscillator. If 



18 



S.A. Chorosavin 



something had reflected the waves, more precise, if the emitted waves had returned 
to the oscillator, we could wait for that the these waves would stimulate an increase 
of the amplitude of the oscillator. But if no wave returns to the oscilatar, -it is 
exactly our case-, the damping dominates and the amplitude decreases This 
phenomenon is said to be a braking radiation or damping radiation. 



5 as we have seen, it decreases to the zero value 
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1.4 Particular Cases. 2. Resonance 

We begin to analyse the inhomogeneous equations with the most simple case, where 
there is an only right-hand incident wave 

uoi(t,0, x) = uo(t, 0,x) = c + (x + t) 

and where one can easy compute the solution: don't deviate from standards, let 

Uoi(i, 0, x) — uo(t, 0, x) = c + (x + t) := Asa\{k(x + tj) 

i.e., 



\dt 2 



(— + 2 7 ^ + ft 2 ) (Q(t) - A sin(fct)) = 2 7 ^|^ - 2 7 A/c cos(fci) 



According to the standard reading of the universiy course to the theory of ordinary 
differential equations, if 

(jjp + 2 "f-Q t + Q2 )y(t) = const + Q 2 Asm(kt) 



/ d 2 d \ 

+ 2 7 — + ft 2 J Y(t) = const - 2 1 Ak cos(fci) 



then a partial solution to these equations is given by: 

»<" = (=i? + Sh <w? (( ' ta + tf) sin(t " ~ 2 ^ cos < t ") + ^ 



On the surface, there is no essential deflecting from the situation in the university 
course of ordinary differential equation, but it is only on the surface. 



dq(t) 

at 



the forms 

t=o 



Through the presence of the term 

q(t) = A s sm(kt) + A c cos(fct) + Const 



need not contain a particular solution to the former equation. Although, we have, 
of course, 

fl 2 A 

q(t) = (_ fc 2 + ^2)2 + (2 7 fc)2 (( _fc2 + fi2 ) sin ( fc *) - 2 l k cos(fci)) + constx 



+ 2 7 dq(t) 



t=0 



n 2 at 

+ e~ 7 * (consti cos(ft 7 i) + consti sin(ft 7 i)) 
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we can write this relation as follows: 

?(*) 



^(-fc 2 + n 2 ) 2 + (2 7 /c) 2 



sin(fct + 4>k) + consti 



t=o 



n 2 at 

+ e~ Jt (consti cos(fi 7 t) + const 2 sin(f2 7 i)^ 



where (f>k is such that 



cos (p k 



-k 2 + n 2 



v / (-k 2 + n 2 ) 2 + {2 1 k) 2 



, sin 4> k = - 



2 7 fc 



sj{-k 2 + n 2 ) 2 + (2 7 fc) 2 



As for Q(t) — Asin(fci), on referring to the equation to Q(t) — As'm(kt) we see, 



dg(t) 



dt 



is there exterior and docs not hinder for 



t=o 



on the contrary, that the term 2 7 

immediate displaying a particular solution to the lonely equation to Q(t) — As'm(kt) 
itself. It is this: 



Q(t) - Asin(fci) 



-2jAk 



(2 7 fcsin(fct) + (-fc 2 + Q 2 ) cos(fci)) + 



2 7 dq(t) 



{-k 2 + n 2 ) 2 + {2 1 k) 2 

Nevertheless, this form is improper, because of the initial relations 



n 2 dt 



(Q(0)-c+(0)) =0 
which in our case read 

Q(0) = 0, 



d(Q(t)-c+(t) 



dt 



=o = 2 7 ( g (0)-c + (0)) 



d(Q{t) - Asin(fci)) 



dt 



t=o 



279(0) 



Once again, we need to take into account the terms containing e 7 *(' ' ' ) : 

-2jAk, f*..,.^,^ , , ,.2 , r^_„,^ , 2 7 9g(i) 



Q(i) - Asin(fci) 



(-fc 2 + ft 2 ) 2 + (2 7 fc) 2 
+ e~ 7 * (consti cos(r2 7 t) + consti sin(f2 7 i)^ 

: COs(fci + 



-2 7 Afc 



(2 7 fcsin(fct) + (-fc 2 + f7 2 )cos(fci)) + ^ 
27 0?(t) 



and 



^(-k 2 + n 2 ) 2 + (2 7 fc) 2 ~" v "" 1 1 fi 2 at 

+ e -7 * ^consi 3 cos(ri 7 t) + consti sin(f2 7 t)^ 
Let us restrict ourselves to the case where yt » 1. In this case e -7 *^- "j«0 



?(*) 



ft 2 A 



^(-fc 2 + tt 2 ) 2 + (2 7 fc) 2 
27 dq(t) 



sin(fci + 4>k) + consti 



+ 



n 2 at 



t=o 



Q(i) - Asin(fci) 



2 " ' U (2 7 fcsin(fct) + (-fc 2 + f! 2 ) cos(fct)) + 27 9<?( * } 



(-fc 2 + n 2 ) 2 + (2 7 /c) 2 
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If k 2 — 0? — 27 2 , □ then the amplitude of the harmonic part of q(t) , i.e., the value 
of the quantity 

\n 2 A\ 



becomes maximal @: This phenomenon is said to be a resonance 



^/(-P + 2 ) 2 + (2 7 fc) 2 



dq(t) 



Next, if k = f2 , u and if, in addition, — ^ 
Q(t) - Asin(kt) t 



t=o 
2-yAk 



{2 1 k) 2 
- Asin(fci) 



0, then 
(27fcsin(fct)) 



Let us now recall, that 



u(t, x) = 



Q(t - \x\) - c+(t ~ \x\ 





if <t- \x\ 
if t-\x\ <0<t 
if t < <t + \x\ 



and hence, if fc 2 = Q 2 , ^f- 



Q(t+\x\) + c+(t+\x\) , ift+|x|<0 
0, jt » 1, then 



t=o 



-Asm(k(t - \x\)) , \iQ<t-\x\ 



u(t,x) « 

In particular, 

u(t, x) 



if t-\x\ <0<t 
if t < <t + \x\ 



+ c + {x + t) 



+ Asin(k(x + 1)) 



-Asm(k(t+\x\)) , ift+|a:|<0 





2A sin(kx)cos(kt) 
Asm(k(x + 1)) 



We see, in the case of 



fc 2 = ^ 2 ,^ 
at 



if < t- \x\,x < 0,7< >> 1 
if 0<t - \x\,x > 0,7< >> 1 
if t- \x\ < 0<t 



0, 7* » 1 



the incident wave is essentially completely reflected by the oscillator!!! This phe- 
nomenon can also be referred to as a kind of resonance. 



6 fix on: fc 2 = D. 2 - 27 2 , not fc 2 = Q? !! 
7 up to assumption made just now 
8 fix on: fc 2 = U' 2 , not k 2 = Q 2 - 2 7 2 !! 
9 hence, Q(t) at 
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1.5 Comments. Discussion. 



In this section, we have considered two models of interaction and obtained for 
oscillators the corresponding equations: 



d 2 q(t) 

dt 2 



-2 7 ^ - n 2 q(t) + 27^| t=Q + tfu a (t, 0,0) + f Q (t) + 2jJ* f (r)d7 



d 2 go(t) 

dt 2 



= -n 2 q (t) + 2 l2ll [ gb(r)dT + 71^01 (t,0,0) + / (t) 
Jo 



These equations are not literally ordinary differential equations , ordinary dif- 
ferential equations, because of 'unordinary' terms 



+2 7 



dq(t) 



dt 



t=0 



and, resp., 27172 



/'*>(• 

Jo 



r)dr 



We can reduce these equations to ordinary differential equations, but the latter 
turn out to be third order equations: 



d 3 g(t) 
Ot 3 



-2 7 



d 2 q(t) 

dt 2 



dt at 



! Uo (M,0) + ./o(t) + 2 7 f /o(r)rfr) 



d 3 g (t) 
dt 3 



dt 



+ 71 



(27250 (t) 



+ 



du (t, 0,0) \ , 0/o (t) 



at 



at 



We concentrate now on the system's behaviour at large i . In order to estimate 
the asymptotic behaviour of q(t) and qo{t), as t — > +00, let us handle with the 
characteristic polynomials. They are: 



a 3 + 2 7 a 2 + n 2 \ 

The first polynomial has the roots 

Ai = 0, A 2 = -7+ Vt 5 " 



Aq + fi 2 A - 2 7l72 



-7 



and we see that first root is zero and two other roots have non-positive real parts. 
Hence, except for the unique case, the case of 



+2 7 



Oq(t) 



dt 



t=0 



we have exponentially decaying q(t) as t — > +00. Whatever the case, q(t) is bounded 
as t — > +00. Recall, if /o and uq both are identically zero, then 



q(t) -> 2 7 



at 



t=0 



as t — > +00 . 



However, the limit function 



qoo(t) = 2 7 



&?(*) 



at 



t=0 



is NO solution to 



dt 2 



dt 



dt 



t=0 
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every time that 

2 7 ^ 



?0. 

t=o 



m 

As for the second polynomial, and the second equation, the situation is more 
dramatic, much more. The polynomial 

xl + n 2 \ - 2 7l72 

has one pure real root A i and two complex-conjugated ones: A 2 , A 3 , A 3 = A 2 • 
Since 

Aoi + A02 + A03 = , A01A02A03 = 27172 

we have 

A i + 2i?eA 02 = , A i | A 02 1 2 = 2 7i72 

and hence 

A O i>0 (!!!) , J ReAo2 = iieAo3<0. 

Thus, we observe an EXPONENTIAL GROWTH of the oscillator amplitude, as 
t — > +00, the fact used to perplexing physicist's mind. We cannot here hope we 
have simply confused the 'time direction'. If we had, we would have two roots 
with strictly positive real parts! We defer the more detailed discussion on this 
subject and notice only, that a related phenomenon is known in electrodynimics, 
see Abraham-Lorentz-Dirac equations. 
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A little more general Model with Finite Rank 
Interaction. The Case of non-local Interaction 
and arbitrary Interaction Forces. 



We will now. slightly change the model. We will not suppose that 8 shall be Dirac's 
5-function E3and B be Jp- on the whole line. At the beginning of this section, we will 
consider an abstract d'Alambert-like equation and recall the associated d'Alembert- 
like solution formulae. Then we introduce an abstract analogue of the interaction 
discussed in the previous section. Then we discuss some partial cases, imitating the 
reasons of the same section, and finally indicate how the abstract formulae correlate 
with the formulae obtained in the previous section. 
Now then. When solving the equation 



dt 2 



denote 



Ml 



= Bu + f 

du 
~dt' 



Write the equation to be solved as 

:(:m 



d_ 

dt' 



Ul 

Bu + f 



1 
B 



u 

Ul 



Let 



V, 



denote the propagator for 



dV 



V s 




' 






Vi j 


) = 


v B 


I) 


C) 



Then 



( «(*) \ 
V m(t) ) 



V, 





fir) 



Apply usual matrix -form representation 

Vt,a = ( 



Vu{t,s) Visit, s) 
V 2 i(t,s) V 22 (t,s) 



Then obtain 



u(t) = Vii(t,s)u(s)+Vi 2 (t,s 



/ dujt) 
\ dt 



Vi 2 (t,T)f(T)dT 



u (t) :=Vii(t,s)u(s)+Vi 2 (t,s)(^^- 



u(t) = u (t,s) 



V 12 {t,T)f(T)dT 



I will now take 



/(*) = S a ,t, Xo U (t,Q) + fi(t) 
Q(t) = < l\u(t) > 



10 In order to no confusion can ocurr we replace the symbol 8 by <5 Q 
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where: 

1) Uo(t, Q) is such that its values must be numbers, and 

2) I is a linear functional. 
In addition I will take 

d 2 q 

In the previous section 5 a ,t,x was a function of x and it did not depend on t 

S a , t (x) = S(x - x ) ; 
As for I, it was such that for any function F of x one had set 

< l\F>=F{x ) 

In addition we had there taken 

U (t,Q) = 4 7c (0(t) - g(t)) , 7c = 7 c, 

F p/ (g, Q, t)/M = -ft 2 (g(t) - Q(i)) + / (t) . 
Now we are mimicking, imitating arguments of the previous section: 

u(t)=u (i)+ / ^12 (t,T)f(T)dT 
J s 

u{t)=u (t) + V 12 (t,T)(~ 5 a , T , Xo Uo(T,Q) + hir^dr 

U(t) = U (t) + J Vvz{t,T)h{T)dT - j V 12 (t,T)(6 a , T>Xo Uo{T,Q)yT 

u i(t) :=«o(*)+ / Vi 2 (t,r)/i(r)dT 

= «oi(*) - ^ Vl2(£,T)(£ a , r , X0 £^(T,Q))dT 
< i|«(t) > = < i|«oi(t) >- <l \f V 12{t,T)(S a , T , X0 U {T,Q)^dT > 

We will suppose 

< l\ J* V 12 (t,T)(6 atTtX0 U (T,Q)yT >= J* < l\V 12 (t,T)(6 atT , X0 U (T,Qj) > dT 
Hence 

< l\u(t) >=< l\u i{t) > - I <l\V 12 {t,T)S a . T . Xa > U (T,Q)dT 

J s 

Recall 

Q{t) :=< l\u(t) > 

and put 

K (t,T) :=<l\Vi 2 (t,T)6 a , T >, Q (t) :=<l\u (t) >, 
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Qoi(t) :=< J|uoi(t) >=< i|uo(t) > + f < l\V 12 {t,T)h(T) > dr 

J s 

Then obtain 

Q(t) = Q (t)- J K (t,T)U (r,Q)dT 
Let us now define an operator Kq\ 

(k h)(t) := f K (t,T)h(r)dT 

J s 

Then the recent relation becomes as following: 

Q = Qoi-Ko(u (-,q)) 

If we take into account that 

Mg£=F pf ( q ,Q,t) 
then the equations of motion become 

Q = Qoi -#o(l>b (•,<?)) 

The character of this relation array is too general. Little can be said about 
its properties without any specifying. So, let us restrict ourselves: First, we will 
assume 

U (t,Q) = 4 7c (Q(t) -<?(*)) 
F pf (q,Q,t) = -K(q-Q)+F (t)=M(-n 2 (q-Q) + f (t^ 

and denote 




Then the equations of motion become 

Q = Qoi - 4 7c ^o(Q - q) 

q = -n 2 (q-Q) + f 

We infer: 
1st step, 

Q + 4^ C K Q = Qoi + 4-f c K q 

q = -n 2 ( q -Q) + f 

2nd step, 

(l + 4 7c ^ ) Q = Qoi + 4j c K Q q 

q = -n 2 q + n 2 Q + f a 
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3rd step, 

(l + 4 7c ^ ) Q = Qoi + ^lcK Q q 
(l + 4 7c K )g - (l + 4 7c ^ )(-f! 2 ( z + f! 2 Q + ./o 

4th step, 
(l + 4 7c if ) Q = Qoi + 4 lc K q 

(l + 4 7c X )g = -{l + 4 lc K )n 2 q+ (l + 4 7c X )r! 2 Q+ (l + 4 7c £ ) /< 

= -n 2 (i + 4 7c x ) q + n 2 (i + 4 7c ir ) Q + (i + 4 7c ir ) /o 

replace here ^1 + 4 7c i\ ^(3 
by (Qoi + 4 7c iT 9) 
= -ft 2 (l + 4 7c AT ) q + ft 2 (Qoi + 4 7c A r) <?) + (l + 4 7c AT ) /o 
- -fl 2 q + n 2 Q i + (l + 4 7c ^ ) fa 

We have now seen: 

(i + 4 7c x ) q = -n 2 q + n 2 Q Q1 + + 4 7c ir ) f 

It means 

= -n 2 q(t) + n 2 Q 01 (t) + f (t)+A lc J^ Ko(i,r)(/o(r))dT 
We rewrite this equation as following: 

^ = -4 lc J' Ko(t,T)^^dr-n 2 q + n 2 Qoi(t) + fo(t) + 4 lc J' Ko(t,T)fo(r)dT 
and then, replacing Qoi(i) by < Z|u i(i) > , as 

+ ft 2 < ZKi(t) > +f (t) + 4 7c / K (t,T)f {r)dT 

J s 

Thus, we have obtained an insulated equation for q, and the equation resembles 
ones in the previous section. To continue, we must specify the nature of the abstract 
terms, of u{t). 

If they, u(t), all are elements of an functional space, say simple, if u(t) are 
functions of a 'spatial' variable x, our constructions and formulae appear as follows: 

u{t,x) = u (t,s,x) + J {Vl 2 (t,T)f(T)^(x)dT 

u(t,x)=u {t,s,x) + J J Vi 2 (t,T,x,£)f(T,£)d£dT 
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In this case we have 

f(t,x) := -S (x)Uo(t, Q) + fi(t, x) 
( in the previous section, 5 a j, Xo (x, x ) = S(x — x ) ) Then 

u{t,x) = u (t,s,x) - J V 12 (t,T,x,t)(5 atTtXo (£)U (T,Q) - /i(T,£))d£efrr 
Hence 

< l\u(t) > = < l\u(t, •) > 

= <l\u Q (t,3,-)- J J V 12 (t,T,-,0(S a>T , Xo {OUo{T,Q)-f 1 {T,0)d^dT> 
= <l\uo(t,8,-)>-<l\J J^V 12 (t,T,;t)6 atTtXO {Z)U Q {T,Q)d£dT> 



Is JX 

+ <l\j j V 1 2(t,T,;0f 1 (T,£)d£dT> 
ft 



<l\u (t,S,-)>- <l\ V 12 (t,T,;!;)5 
Js JX 

+ <l\J*f v 12 (t,T,;Qf 1 (r, £)d$dT > 



Hence 

< l\u(t) > = < l\u(t, •) > 

= < 



i.e. 



Of course 



l\u {t,s,-) > -J < l\J V 12 {t,r, ;06 a , T , X0 (Z)dt > U (T,Q)dT 
Q(t) = QoiW- f <l\ [ V 12 {t, t, ;Q5 a , T , X0 (Qdt >U (r,Q)dr 

Js JX 

Q(t) = Qoi(i)-/ K (t,T)U {T,Q)dT 

J s 

K (t,r) :=< l\ [ V 12 (t,T,;06 a ,r >X0 (0<% > 

JX 



IX 

r-l 



Q 01 =<l\ U0 {t,S,-)> + <l\ [ f V 1 2(t,T,;Z)f 1 (T,S)d£dT> 

Js JX 

If we presume B, I to be the same as before, in the previous section, i.e., if B is 
c 2 ^- on the whole line, and < l\F >= F(x ), then 



' V 12 {t,T,X,S)f{T,S)d€dT 
Is JX 

ft 



\ f (/(r, x + c(t- r)) - />, x-c(t- r))) dr 



j j V 12 (t, t, x, OS a ,r, Xo (OU (t, Q)d£dr 



Y c j (^a,r, Xo {x + c(t-T))-6 a ^ Xo (x-c(t-T))jU (T,Q)dT 
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Thus we obtain 

K (t,T) = ^(s a ,T,x (xo +c{t - r)) - S a ^ Xo (x - c(t - t)) 
Here S a _ TtXQ stands for any primitive of <5 QjTjXo , i.e., 

Note once again that 

Sa,r,x (x + C(t - t)) - S a . T . Xa (x - c(t - t)) 

does not depend on whatever primitive of 5 a ,T,x a which one has chosen!!! 
Suppose, S a . t ,x ( x ) is °f the form 

$a,t,x (x) = S a {x - x ) 

Then we obtain 

K (t,r) = ^(5 a ,T,x (xo + c(t - t)) - S a . T . Xo (x - c(t - r))) 

= ±(6 a {c{t-T))-6 a {<T-t)j) 

In this case, we have finally seen, 



j)j2 = -2^ J (5 a (c(t - r)) - L(c(t - i))) ^j^dr - ft 2 g + n 2 Q i(t) 

+ f (t) + 2^ J (5 a (c(t - t)) - ~S a (c(T - t)))/ (r)dr 
= -2jJ* (S a (c(t - r)) - ~6 a (c(r - *))) ^f^r - fi 2 g + ft 2 Q 01 (t) 
+ / (t) + 2 7 ^ (5 a (c(t - r)) - ^(c(r - i)))/o(r)dr 

In any case, 

^ /o(i) + 2^ jf* 2^ (i,r)(^§^ - f (r))dr = -fi 2 (g(t) - Q i(*)) 



dt 2 

or, more generally 



Q(t) = Qoi(*)- / i^o(t,r)(7o(T,0)dT, 
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3 APPENDIX 



3.1 APPENDIX A. Abstract Linear Response Formula 

We recall some abstract linear response formulae. 

First, we recall linear response formula for the first order differential equation. 
We indicate only formula itself, minimum of details. We say only, a deriving of 
this formula is based on an abstract modification of the method of variation of 
constant. As for detailed description, how one produces such a formula, as well as 
for proofs, background and all that, see suitable standard manuals. 

Let {A(t)} t stand for a family of linear operators, so that the homogeneous 
linear equation 



dv(t) 
dt 



= A(t)v{t), v(t) 



^initial i 



(t>s) 



has a unique solution v(t), whatever 'initial time' s and initial data Viniual may be. 
Let 

Vt.s 

denote the propagator, alias evolution operator, to the former system, i.e., 

V(t) = Vt.sVinitial = Vt.sV(s). 

One can show that 



Vt. s , if t > r > s, ( consistency relation ) 



V t . t = V r ,r = V s , s = I . 

The linear response formula reads: The formula 



U(t) = Vt t sUinitial + / Vt, r /(r)dT 
J s 



gives a solution to the inhomogeneous equation 



du(t) 
dt 



= A(t)u + f(t), u{t) 



^(^) — ^initial- 



We will write the response formula as following: 



U{t) = V t>s u( S ) + f V tlT f{T)dT. 
J s 



One can now produce the response formula for the second order differential equa- 
tions. As befor, the response formula accommodates the formula, which describes 
solution to a homogeneous equation, to the case, where one need to solve the 
associated inhomogeneous one. 
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Now then. When solving the equation 



d 2 u(t) . .du(t) „, . . . „, . 



dt 2 



denote 



dt 



ui(t) := 



Write the equation to be solved as 



d_t u(t) \ f 
dt\ J 



du{t) 
dt ' 



m(t) 



K(t)u { t) + B(t)u(t)+f(t) 

/ 1 W U (t) N / N 

V B(t) A"(t) A ui(i) J + I /(t) J 



Let 



denote the propagator to 



Vt 



t,s 



d_( v(t) \_(0 1 w v(t) \ 

J I s(t) A-(t) A «i(t) J 



Then 



Apply usual matrix-form representation 



( 



V n (t,s) V 12 (t,s) 
V 21 {t,s) V 22 {t,s) 



) 



Then obtain finally 



u(t) = Vu(t,s)u(s) + V 12 (t,s){^-\^ s ) + f V 12 (t,T)f(r)dr 



Notice, 



and the formula 



/ Vu(t,s) V 12 (t,s) \_f I \ 
V V 21 (t,s) V 22 (t,s) ) \ I ) 



v(t) = V 11 (t,s)v(s) + V 12 (t,s)( 



dv(t) 



dt 



is a formula which describes solutions to the homogeneous equation. Sometimes, it 
is instructive to notice that 



and hence 



V 21 {t,s) = —V u (t,a), V 22 (t,s) = —V 12 {t,s). 
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3.2 APPENDIX B. 

Response Formula for y(t) + 2^y{t) + Q 2 y{t) = f(t) 
If one applies the response formula to the equation in the title, then one finds: 

V(t) = e-*-> (cosO^ - s) + 7 Sin "^ = S) )y(s) 
l{t -s) smOyft - s) / dy{t) 

fi 7 I at 



e 



7(t _ T) 8 mn^t-T) /(T)dT 



where ft 2 = _ ^ [u] 

Actually, the solution to the joined homogeneous equation 

x(t) + 2 7 ±(t) + n 2 x(t) = 

is this: 

sinfLix , . _^ t smQ~,t 



= e~ 7t (cos» 7 t + 7 7 WcO + e" 7 * 7 x(0) 
x(t) = e - 7t (-fi 7 2 _ 7 2 )(5^) x( o) + e - 7 *(cosfi 7 i- 7 ^^)i(0) 



7 ' v "7 

Since SI 2 = f2 2 — 7 2 , we can besides write 

sinfLix , . sin f2.-»t 



x(t) = e~ 7t (cos» 7 t + 7 7 )a:(0) + e~ 7t 7 x(0) 
i(t) = e - 7 *(-fi 2 )(^^)x(0) + e - 7t (cosO 7 t-7^^)i(0) 



The coefficients of the latter equation for x{t) are constant. Hence 
x(t) = e-^-^ (cosil^t - s) + 7 Sin "^" s) )x( S ) 

+e _ 7(t _ a) sinn^- a ). (a) 

»(t)=e-T(«-)(-tf) Bmi y- a) x( a ) 



sin f2 7 (t — s) 

7 

+e- 7 (*- s ) ( cos tt 7 (t - s) - 7 Sin ^~ S) ) 



Thus we have seen, that the formula for the propagator Vt. s , which associates with 
the homogeneous equation x(t) + 2 r yx{t) + il 2 x(t) = 0, is this: 



Vt.s = 



V n (t,s) V 12 (t,s) 
V 21 (t,s) V-afaa) 



(COS n T (t - S) + 7 Sln y 3) ) e - 7 (t- s) sin0 7 (t- S ) 

e- 7 (*- s )(-fi 2 )^%^ e -7(t-«) (cos0 7 (t - a) - 7 sin y } ) 



11 if = 0, then read s "^, T = t, of course 
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u(t) = V u (t,a)u(a) + V 12 (t,a) 



and replace there the proper terms according to 



V n (t, s) = c-r(*-) (cosO^ - s) + 7 ™Mz£)) 



F 12 (M) = e-^) Sin ^" S) 
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3.3 APPENDIX C. 
Particular Solution to 

y(t) + 2>yy(t) + Vt 2 y{t) = A s sin(H + 0) + A c cos(H + 0) 
We have 

(^2 + 2r Tgl + ^ ( a s Sill ( kt + 4>) + a c COs(kt + 

= (~^ 2 + fi 2 ^ a s sin(fci + 4>) + 2jka s cos(fci + </>) 
- 2jka c sin(fci + </>) + (-fc 2 + Q 2 ^j a c cos(fci + <j>) 

In order to 

/ -k 2 + Q 2 -2 7 fc \ ( a s \ ( A.' 
{ 2 7 fc -k 2 + n 2 ){a c J { A c 

it is sufficient that 

((-fc 2 + n 2 ) 2 + (2 lk f) ^*s\_f -k 2 + n 2 2 7 fc ^ ( a s 

Thus, if 



-2 7 fc -fc 2 + Q 2 J \ A c 

(-fc 2 + n 2 ) 2 + {2 1 kf ± o 



then 

/ a s \ 1 / -fc 2 + Q 2 2 7 fc \ ( A s 



a c J (-k 2 + fl 2 ) 2 + (2 7 fc) 2 \ -2 7 fc -fc 2 + tt 2 J \ A, 
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3.4 APPENDIX D. Helpful formulae 

Derivatives. 1) 

de -7 ' (A cos(ut + ip) + B sin(wt + ip)j 

dt 

= -7e~ 7t (a cos(ut + ip) + B sin(o>t + ipfj 

+ e~ 7 * (^— loA sin(wt + ip) + luB cos(ujt + ipfj 
= e~ 7 * ((-7A + ujB) cos(ut + ip) + {-ujA - jB) sin(wt + ipfj 

2) 

<9e -7 * ^— 7 cos(wt + ip) + lu sin(wt + ipfj 
dt 

3) 

<9e~ 7 * (uj cos(wt + ip) + 7 sm(ujt + ipfj 

dt = 

2a) 



e 7 *(7 2 + uj 2 )cos(ut + ip) 



e" 7 *(7 2 +uj 2 )sm(ojt + (p) 



g e -7(t-T) cos ^(t - T ) + p)+uj sin(w(i - r) + 



<9t 

= - e " 7 (*-% 2 + cj 2 ) cos(u;(t - t) + if) 



3a) 

<9e~ 7 * cos(w(t — t) + </?)+ 7 sin(w(t — r) + ipj*J 
dr 

= e - 7 (*- T )( 7 2 + to 2 ) sin(w(t - r) + 

Other. 

2 sin a sin /? = cos(a — (3) — cos(a + (3) 

U)(t — t) - UJ S (T — to) = (LU+LU s )(t — T)—LJ s (t — T)—LJ s (T — to) 
= (w + 0J s )(t - t) - U> s (t - t Q ) 



u(t - t) +uj s (t - t ) 



= (u> - u s )(t - r) + UJ s (t - t) +u s (t - t ) 

= (UJ - U s )(t - T) + LJ s (t - t ) 



2 sinw(t — t) sino) s (r — to) 

= cos((cj + UJ s )(t - t) - u> s (t - t )) - COS^(o) - LU s )(t - t) + ui s (t - t )) 
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